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Context. Neutrinos heavier than M z /2 ~ 45 GeV are not excluded by particle physics data. Stable neutrinos heavier than this might 
contribute to the cosmic gamma ray background through annihilation in distant galaxies as well as to the dark matter content of the 
universe. 

Aims. We calculate the evolution of the heavy neutrino density in the universe as a function of its mass, M N , and then the subsequent 
gamma ray spectrum from annihilation of distant NN (from < z < 5). 

Methods. The evolution of the heavy neutrino density in the universe is calculated numerically. In order to obtain the enhancement 
due to structure formation in the universe, we approximate the distribution of N to be proportional to that of dark matter in the GallCS 
model. The calculated gamma ray spectrum is compared to the measured EGRET data. 

Results. A conservative exclusion region for the heavy neutrino mass is 100 to 200 GeV, both from EGRET data and our re-evalutation 
of the Kamiokande data. The heavy neutrino contribution to dark matter is found to be at most 15%. 

Key words. Elementary particles - Neutrinos - (Cosmology:) dark matter - Gamma rays: observations 
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1. Introduction 

The motivation for a fourth generation neutrino comes from the 
standard model of particle physics. In fact, there is nothing in the 
standard model stating that there should be exactly three gener- 
ations of leptons (or of quarks for that matter). 

The present limits on the mass of a fourth generation of neu- 
trinos are only conclusive for < Mz/2 w 46 GeV (Yao et al. 
2006, p. 35). This limit is obtained from the measurement of the 
invisible width of the Z°-peak in LEP, which gives the number 
of light neutrino species, as N v = 2.9841 ± 0.0083 (The LEP 
Collaborations 2001). 

In Maltoni et al. (2000), a fourth generation of fermions is 
found to be possible for M N ~ 50 GeV, while heavier fermions 
are shown to be unlikely. However, this constraint is only valid 
when there is a mixing between the generations (Novikov et al. 
2002) and since this is not necessarily true, we will not take it 
for certain. 

In the context of cosmology and astrophysics there are other 
contraints. Light neutrinos, with 55 1 MeV, are relativistic 
when they decouple, whereas heavier neutrinos are not. The light 
neutrinos must have 2 m v < 46 eV in order for Q. v h 2 < 1 to be 
valid (Hannestad 2006b). For the dark matter (DM) content cal- 
culated by Spergel et al. (2003), the bound is £ m v < 12 eV. 
The number of light neutrino species are also constrained to 
N v = 4.2*J - 2 by the cosmic microwave background (CMB), large 
scale structure (LSS), and type la supernova (SNI-a) observa- 
tions at 95% confidence (Hannestad 2006a). 

Neutrinos heavier than about 1 MeV, however, leave thermal 
equilibirum before decoupling and therefore their number den- 
sity drops dramatically, see for example Dolgov & Zeldovich 
(1981). This will be discussed in more detail in Sect. 2. 

The most important astrophysical bound on heavy neutrinos 
comes from Kamiokande (Mori et al. 1992) and this will be con- 
sidered separately in the end. 
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In Fargion et al. (1995), it is found that the mass range 60 < 
< 1 15 GeV is excluded by heavy neutrino annihilation in 
the galactic halo. However, according to Dolgov (2002, p. 57) 
this constraint is based on an exaggerated value of the density 
enhancement in our galaxy. 

Other works constraining the heavy neutrino mass include 
Fargion et al. (1998, 1999) and Belotsky et al. (2004). There 
has also been a study of the gamma ray spectrum of dark matter 
(DM) in general (Ando et al. 2007). 

For an exhaustive review of modern neutrino cosmology, 
including current constraints on heavy neutrinos, see Dolgov 
(2002). It is concluded that there are no convincing limits on 
neutrinos in the mass range 50 < < 1000 GeV. A review of 
some cosmological implications of neutrino masses and mixing 
angles can be found in Kainulainen & Olive (2003). 

In this paper we consider a stable fourth generation heavy 
neatrino with mass Z 50 GeV possessing the standard weak 
interaction. We assume that other particles of a fourth generation 
are heavier and thus do not influence the calculations. 

We assume a ACDM universe with Cl tot = Q,„ + £2a = 1, 
where Q,„ = Q b + Q. DM = 0.135//i 2 , Q. h = 0.0226 / 'h 2 and h = 
0.71 (Spergel et al. 2003), using WMAP data in combination 
with other CMB datasets and large-scale structure observations 
(2dFGRS + Lyman a). 

Throughout the article we use natural units, such that the 
speed of light, Planck's reduced constant and Boltzmann's con- 
stant equal unity, c — h = ks = 1. 

If heavy neutrinos (Mpt > 50 GeV) exist, they were created 
in the early universe. They were in thermal equilibrium in the 
early stages of the hot big bang, but froze out relatively early. 
After freeze-out, the annihilation of AW continued at an ever 
decreasing rate until today. Since those photons that were pro- 
duced before the decoupling of photons are lost in the CMB, 
only the subsequent AW annihilations contribute to the photon 
background as measured on earth. 
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The intensity of the photons from AW-annihilation is af- 
fected by the number density of heavy neutrinos, n^, whose 
mean density decreases as R 3 , where R is the expansion fac- 
tor of the universe. However, in structures such as galaxies the 
mean density will not change dramatically, and since the number 
of such structures are growing with time, this will compensate 
for the lower mean density. Note that the photons are also red- 
shifted with a factor R due to their passage through space-time. 
This also means that the closer annihilations will give photons 
with higher energy than the farther ones. 

2. Evolution of neutrino density 

Let us recapitulate the results of Dolgov & Zeldovich (1981). 

The cosmic evolution of the number density, tlx, of a particle 
X can, in general, be written as 



h x = -n\ (cry) - 3H(t)n x + i//(t), 



where (crv) is the thermally averaged product of the mean ve- 
locity and total annihilation cross section for the particle, and 
H(f) = R/R is the Hubble constant. The term -3H(t)rtx repre- 
sents the expansion of the universe, and the production term is 
i//(t) = n Xeq (crv), where nx eq is the equilibrium concentration of 
particle X. 

If we write r x = n x /n y , Eq. (1) can be expressed as 



r x 



(crv) n y (r x 



where 



fXeq 



(2ny 



' Xeq 
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if 6 = T/m x > 1 

ife< l. 



(2) 



(3) 



Here f(3) ~ 1.2020569 is the Riemann zeta function, T is the 
temperature, nix is the mass of particle X and g s is the number of 
spin states. For photons and electrons, g s = 2, while formassless 

left-handed neutrinos, ? t = 1. For reference, , (2 "1, , » j. 

The value of the relative equilibrium concentration, rx eq , is 
derived from 
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where the expressions for n y and n eq were taken from Dolgov 
(2002, Eq. 30). 

According to Dolgov & Zeldovich (1981, Eq 2.9), freeze- 
out (equilibrium destruction) occurs when the rate of change 
of the equilibrium concentration due to the temperature de- 
crease is higher than the reaction rates, which means that 
2 (crv) n y rxeqtT I m > 1. Until freeze-out, the relative parti- 
cle density follows the equilibrium density closely: rj X ~ r Xeq - 
Hence, the relative density at the moment of freeze-out is 



r fx = (2 (crv) n y t f f ) 



rxeq, 



(5) 



where tf and Of - Tf/trtx are the time and relative temperature 
at freeze-out. 

As the temperature decreases, the production term r Xeq will 
drop exponentionally, such that the relic concentration of X will 
be more or less independent of rx eq - With this approximation 
( r Xeq = 0), Eq. 2 can be solved for t — > oo: 
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where we have used tT 2 » 3.677 x 10 18 / yfgl (Dolgov 2002, 
Eq. 37), with g,(T f ) from Kolb & Turner (1990, p. 65) be- 
ing the number of relativistic species in thermal contact with 
the photons. Furthermore, n y (to) = 2T^(3)/tt 2 m 0.247"g is 
the photon density today, where the photon temperature today 
is T = 2.725 K (Mather et al. 1999). According to the stan- 
dard model of particle physics, g* = 106.75 for T > 100 GeV 
(g* w 200 for supersymmetry models at yet higher tempera- 
tures). If we assume that Of «: 1 (which we will later show to be 
reasonable), we obtain r^x ~ rfxOf, which differs by a factor two 
from the result of Dolgov & Zeldovich (1981, Eq. 2.11). This is 
natural if they consider the density of n N+ fj since our rox is valid 
for N and N separately. 

In order to take into account the increase in temperature due 
to entropy conservation after freeze-out of particle X, we must 
take 



43/11 

rOX 7^7T»y( f 0) 
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(7) 



(In fact gj f /2 should be written g^ f ■ g x J f 2 but for T f > 0.1 GeV, 

g*Sf = g*f.) 

We now turn to the case of heavy neutrinos. Since we wish 
to avoid the lenghty calculations of the cross sections of heavy 
neutrinos (Enqvist et al. 1989), we use Fargion et al. (1995, 
Fig. 1 and Eq. 4) and solve for (crv). We assume that they use 
g* - g*(Tf) x g»(MA(/30), but the exact value does not change 
the result in any significant way. The resulting (crv) is presented 
in Fig. 1. The cross section drops from ~ 45 GeV, where the 
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Fig. 1. The cross section times the velocity (in m 3 /s) of heavy 
neutrino annihilation NN as a function of their mass (in GeV) at 
freeze-out, T = Tf. 

Z° resonance peaks until the annihilation channel starts 

to dominate at M N > 100 GeV. 

According to Fargion et al. (1995), the cross sections of 
heavy neutrinos can be estimated using the annihilation chan- 
nels 



NN 
NN 



ff 



(8) 
(9) 



There are several other possible annihilation channels for NN — > 
W + W-, like NN -> LL, H°H°, Z°Z° -> W + W~ and also in- 
terference between L and Z°, as well as between L and H°. 
However, in the limit s — > 4M^, which is valid for cosmolog- 



2 (crv) riytf (1 + Of) 2 (crv)/ « r / , - T 'j. (1 + Of) j ca j heavy neutrinos, the dominant channel is through s-channel 
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NN -> Z° (Enqvist et al._1989, p. 656). Furthermore, the other 
annihilation products, AW — > H°H°, Z°Z°, are suppressed with 
respect to -production (Enqvist et al. 1989, p. 651, 656). 

Hence, the above estimation of the (crv) should be fairly accu- 
rate. If anything, it is slightly underestimated. 

Using Eqs. 5 and 3, we can solve for Tj = Gj-M. The result is 
presented in Fig. 2. Note that although it looks like a straight line, 
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Fig. 2. The freeze-out temperature (in GeV a; 1.16 x 10° K) of 
heavy neutrinos as a function of their mass (in GeV). 



it really is slightly curved. We notice that Tf/M^ ~ 1 /30, which 
shows our assumption » 7/ to be valid. This is also in 
agreement with previous results, see e.g. Kolb & Turner (1990), 
where a value of Tf/M^ ~ 1 /20 is quoted. 

We now return to Eq. 7 and apply it to the case of a heavy 
neutrino. We plot the resulting relative relic neutrino density as a 
function of the mass in Fig. 3 using Q/v = 2Mpj ■ «iv(7b)/Pc> 
where p c m 9.47 x 10~ 27 kg/m 3 is the critical density of the uni- 
verse. The resulting heavy neutrino density is very similar to the 
one obtained by Fargion et al. (1995, Fig. 1). The numerical sim- 
ulation also shown in the figure will be the subject of the next 
section. 
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Fig. 3. The relic relative density of heavy neutrinos as a function 
of their mass (in GeV). 



3. Numerical simulation of the neutrino density 

For comparison, we evaluate the evolution of the heavy neutrino 
density numerically. Eq. 1 can be rewritten in terms of the tem- 
perature, T: 

pH(T)n(T) + (crv(T)} (n(T) 2 - n eq {T?)\ , (10) 

where 



n eq (T) = r eq n y = 2T 



(Mn_\ V2 



-U S IT 



{T < M N ) 



\2nT! 

and the relation between time and temperature is given by 

dt_ _ -1 / 1 dg tS /dT 
dT " H(T)[f + 3g tS 



(11) 



(12) 



Here the Hubble constant is H(T) = Hq -\JQ.(T), where the total 
relative energy density of the universe is 



Q.(T) = n R (T) ■ R 4 + Q.m ■ R~ 3 + Q. k 



/T 2 + Q A . 



The curvature term O/. 

gJT) 



and the radiation density is 



8*(T ) 



(13) 



(14) 



due to the reheating as particles freeze out. The reheating also 
means that R = g^To/T (Kolb & Turner 1990, p. 68). The 
number of relativistic species still in thermal contact with the 
photons, g*s (T), is given in Coleman & Roos (2003, Fig. 1). For 
the critical region 0.15 < T < 0.30 GeV their Eqs. 8-9 have 
been used to calculate dg*$ /dT. This updated value of g m $ (T) is 
needed to evaluate dg*$ I dT properly. 

Using a fifth-order Runge-Kutta method with adaptive step- 
size control, taken from Numerical Recipes (Press et al. 1992, 
Ch. 16.2), we solve for n(T) in Eq. 10 using the initial con- 
dition n, = n eq {Ti = Mai/ 15), which is well within the re- 
gion of thermal equilibrium for the heavy neutrinos. The result- 
ing relative relic neutrino density is presented in Fig. 3, where 
Qjv = 2Mn ■ rtN(To)/p c as before. We notice that the peak of the 
curve is Q.n(Mn = 140 GeV) m 0.04, which would then account 
for ~15% of the dark matter content of the universe. 

For comparison, we plot the number density of heavy neutri- 
nos (in m~ 3 ) as a function of T for masses 50, 70, 90, 150, 500 
and 1000 GeV in Fig. 4. As we can see, the transition between 
thermal equilibirum density and completely decoupled neutrino 
density is not sharp. This is one of the reasons for the differ- 
ence between the analytical and the numerical relative density in 
Fig. 3. Another reason for the difference is the inclusion of the 
change in g„s in the evaluation of dt/dT. The evolution of g„s (T) 
is the cause of the small "knee" in Fig. 4 seen at T ~ 0.2 GeV 
(the reheating from the quark-hadron transition). Furthermore, 
when electrons fall out of thermal equilibirum at T ~ 1 MeV 
there is another small knee, reducing again the heavy neutrino 
density somewhat. 



4. Dark matter simulations 

In Sect. 3, we calculated the mean density of neutrinos in the 
universe as a function of redshift and the mass of the heavy neu- 
trinos. However, the neutrino annihilation rate, and thus the in- 
tensity from their gamma spectrum, is proportional to the square 
of the neutrino density. This means that inhomogeneities in the 
universe will tend to enhance the gamma ray signal. 
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Fig. 4. The number density of heavy neutrinos (in m 3 ) as a 
function of T for masses 50, 70, 90, 150, 500 and 1000 GeV (in- 
creasing from left to right in the upper right corner). The dashed 
vertical lines represent the calculated value of Tj in Fig. 2. 
Below T = 0.01 GeV, the curves evolve as (T/Tq) 3 ■ g*$. 



In this section we describe how we calculate the inhomo- 
geneities as a function of space and time, assuming only gravita- 
tional interaction between the dark matter consisting of heavy 
neutrinos and other DM particles. The clumping factor (also 
known as the boost factor) can then be used to calculate the ac- 
tual intensity 



dl_ _ c ,qMo 
dz dz ' 



(15) 



where dl^/dz is the intensity contribution from redshift slice dz 
for a homogeneous universe and C(z) is the enhancement due to 
the clumping at redshift z- 

The clumping factor has been calculated in different settings 
before, ranging from Berezinsky et al. (2006) for local cluster- 
ing giving a clumping factor of ~ 5 to Diemand et al. (2005) for 
mini-halos giving a clumping factor of two orders of magnitude. 
For a discussion about the accuracy of approximating the en- 
hancement with a single clumping parameter, see Lavalle et al. 
(2006), though they focus on antiprotons. 

The spatial and temporal distribution of DM in the universe 
is calculated with the GallCS program. The cosmological N- 
body simulation that we are referring to throughout this paper is 
done with the parallel tree-code developed by Ninin (1999). The 
initial mass power spectrum is taken to be a scale-free (n s - 1) 
one, evolved as predicted by Bardeen et al. (1986) and normal- 
ized to the present-day abundance of rich clusters with cr 8 = 0.88 
(Eke et al. 1996). The DM density field was calculated from 
z = 35.59 to z = 0, giving 100 "snapshots", spaced logarith- 
mically in the expansion factor. 

The basic principle of the simulations is to distribute a num- 
ber of DM particles N 1 with mass M DM in a box of size L 3 . Then, 
as time passes, the particles interact gravitationally, clumping to- 
gether and forming structures. When there are at least 20 parti- 
cles together, it is considered to be a DM halo. It is supposed 
to be no other forces present than gravitation, and the boundary 
conditions are assumed to be periodic. 

In the GallCS simulations the side of the box used was L = 
100/z -1 Mpc, and the number of particles was set to 256 3 , which 
implies a particle mass of ~ 5.51 x 1Q 9 h~ l M Q . Furthermore, for 
the simulation of DM, the cosmological parameters were set to 
D.A = 2/3, Q. m - 1/3 and h = 2/3. The simulations of the 
DM were done before the results from WMAP were published, 



which explains the difference between these parameters and the 
values used elsewhere in this paper, as stated in the introduc- 
tion. Nevertheless, the difference is only a couple of percent and 
should not seriously alter the results. 

Between the initial halo formation at z ~ 11 and the cur- 
rent epoch in the universe, there are 72 snapshots. In each snap- 
shot a friend-of-friend algorithm was used to identify virialized 
groups of at least 20 DM particles. For high resolutions, it is 
clear that the mass resolution is insufficient. Fortunately, the first 
20-particle DM clump appears at z = 11.2, while the bulk of the 
clumping comes from z < 5, where the lack of resolution is no 
longer a problem. 

In order to make a correct large-scale prediction of the distri- 
bution of the DM, the size of the box would have to be of Hubble 
size, i.e., ~ 3000/i _1 Mpc. However, for a given simulation time, 
increasing the size of the box and maintaining the same number 
of particles would mean that we lose in mass resolution, which is 
not acceptable if we want to reproduce a fairly realistic scenario 
for the evolution of the universe. 

We will make the approximation that our single box, at dif- 
ferent time-steps, can represent the line of sight, and since we 
are only interested in the general properties of the dark matter 
clumping, this approximation should be acceptable. 

4. 1 . Validity of simulation 

GallCS is a hybrid model for hierarchical galaxy formation, 
combining the outputs of large cosmological N-body simula- 
tions with simple, semi-analytic recipes to describe the fate of 
the baryons within DM halos. The simulations produce a de- 
tailed merging tree for the DM halos, including complete knowl- 
edge of the statistical properties arising from the gravitational 
forces. 

The distribution of galaxies resulting from this GallCS sim- 
ulation has been compared with the 2dS (Colless et al. 2001) and 
the Sloan Digital Sky Survey (Szapudi et al. 2001) and found to 
be realistic on the angular scales of 3' < 6 < 30', see Blaizot 
et al. (2006). The discrepancy in the spatial correlation function 
for other values of 9 can be explained by the limits of the nu- 
merical simulation. Obviously, any information on scales larger 
than the size of the box (~ 45') is not reliable. The model has 
also proven to give sensible results for Lyman break galaxies at 
z = 3 (Blaizot et al. 2004). It is also possible to model active 
galactic nuclei (Cattaneo et al. 2005). 

Since it is possible to reproduce reasonable correlations from 
semi-analytic modelling of galaxy formation within this simula- 
tion at z = - 3, we now attempt to do so also for somewhat 
higher redshifts. 

4.2. Clumping of dark matter 

We proceed to calculate the clumping factor C(z). The inhomo- 
geneities of the DM distribution can be calculated using the rel- 
ative clumping of dark matter halos: p,- = pilp mea „, where p mean 
is the mean density of the dark matter in the universe and p, is 
the mean density of DM halo i. 

As matter contracts, the density increases, but since the 
gamma ray emitting volume also decreases, the net effect is a 
linear enhancement from the quadratic dependence on the den- 
sity. This means that the DM halos will emit as: 



-halos 



aloi 



(16) 
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where Iq is the intensity for a homogeneous universe and the 
summation is done over all DM halos and thus m, = m na i os . 
The factor Chaio accounts for the modification from the form 
and properties of the halo itself. A simple conic DM distribu- 
tion would give Chaio = 1.6. The more realistic distribution 
p(r) = po • [(1 + + r 2 )]~ l , where r is the radial coordinate 
relative to the halo radius, gives Chaio = l-l- However, the ra- 
diation from within the denser part of the halo will also be sub- 
ject to more absorption, and so for the sake of simplicity we use 
Chaio - 1 ■ We notice that the average relative density over all the 
halos in the simulation is fairly constant, {pi) ~ 70 for z < 5. 

Simultaneously, the DM background (the DM particles that 
are not in halos) will decrease, both in density by a factor 
(rritot - mhaios)/m, ot and because of their decreasing fraction of 
the total mass in the box m tot : 



L DM -background 

T 



Wltot m halos 



m„ 



(17) 



This means that the total clumping factor is 



hialos ^ I DM -background 



V., m,pi | ^ m tot - mhaioA 



m to , 



where the first term starts as unity whereafter it decreases and 
quickly becomes negligeable with respect to the second term, 
which starts at zero, but then rapidly increases. The total clump- 
ing is plotted in Fig. 5 along with the competing (n;v/m~ 3 ) 2 ef- 
fect, as well as the product, all as a function of the redshift z- The 
number density of heavy neutrinos in the figure is taken for the 
mass Mn = 150 GeV. We notice that the clumping enhancement 
remains ~ 30 for z < 1 and that the clumping is ~ 1 for z > 5. 
This is mainly due to the proportion of mass within the halos 
compared to the total DM mass. The clumping enhancement lies 
between the two extreme values by Berezinsky et al. (2006) and 
Diemand et al. (2005) quoted above. 

In fact, the clumping factor can be even higher if other halo 
shapes are assumed with smaller radii (Ullio et al. 2002). The 
densities in the halos considered in the present work have been 
evaluated at the virial radius. 

We also point out that before the reionization, at z > 5, 
there is absorption from neutral hydrogen in the interstellar 
medium (ISM), also known as the Gunn-Petersen effect (Gunn 
& Peterson 1965). This means that photons from higher redshifts 
will be highly attenuated. For z = 5.3, the emission drops by 
roughly a factor of 10, and for z ~ 6 the opacity is T e jf > 20 
(Becker et al. 2001). Hence, any gamma ray signal prior to this 
epoch would have been absorbed. 

5. Photon distribution from A'/V-collisions 

In order to evaluate the photon spectrum from AA-collisions we 
use PYTHIA version 6.410 (Sjostrand et al. 2006). According to 
Enqvist et al. (1989, Eq. 13) the centre of mass energy squared 
is E 2 CM = AM 2 N + 6M N Tf and Tj « M N /30 as estimated above. 

We generate 100,000 AW events for each mass = 
50, 60, 1000 GeV and calculate the photon spectrum and 
mean photon multiplicity and energy. We assume that NN colli- 
sions at these energies and masses can be approximated by v T v T 
collisions at the same E 2 CM . This is obviously not equivalent, but 
AW cannot be directly simulated in PYTHIA. Nevertheless, with 
the approximations used in calculating (cry), the only difference 
between v T v T and NN collisions (except in the cross section) is 
the f-channel production of through r. However, since 
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Fig. 5. The clumping factor (C, dotted line) compared to the 
competing effect of the decreasing heavy neutrino number den- 
sity squared (n 2 N , dashed line) for M# = 150 GeV and the prod- 
uct of the two (solid line). Different neutrino masses scale as in 
Fig. 4. 



the heavy neutrinos are non-relativistic when they collide, the 
two Ws will be produced back-to-back, which means that the 
inclusion of the f-channel is unimportant. 

In order to verify this, we study the difference in the pho- 
ton spectrum for W decay at and 90 degrees, and despite 
an increasing difference between the two cases, even at M N = 
1000 GeV, the difference is not strong enough to change our con- 
clusions. 

The resulting photon distribution is presented in Fig. 6. We 
note that the photon energies peak at £cm/2, which is natu- 
ral since the decaying particles can each have at most half of 
the centre of mass energy. The curves continue to increase as 
oc E~ l as E decreases further. Note that the noise in the curves 
for lower E is due to lacking statistics for these rare events, but 
it does not affect the outcome of the calculations. We also calcu- 
late the mean photon energy and find it to be E y m 0. 21 Ecm for 
all masses. The curve is normalized such that the integral over 
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Fig. 6. The relative energy distributions of photons from 
collisions for heavy neutrino masses M N = 50, 70, 90, 150, 
1000 GeV. Ecm = 2M^ is the centre of mass energy. 
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is unity. The average number of photons, N y , produced for 
an AA-collision is shown in Fig. 7. The sharp rise in the curve 
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at Mm ~ 100 GeV is due to the jets from the emerging W + W 
production. 
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Fig. 7. The average number of photons produced for an NN- 
collision as a function of heavy neutrino mass M N in GeV. 



u 0.01 




£ [GeV] 



1000 



Fig. 8. Cosmic gamma radiation from photons produced in NN- 
collisions as a function of photon energy for neutrino masses 
M N = 50, 70, 100, 140, 200, 500, 1000 GeV. The dotted line rep- 
resents M N = 50 GeV and the dot-dashed M N = 1 TeV The 
solid lines are the masses in between. The circles represent data 
from EGRET (Sreekumar et al. 1998), with error bar, as derived 
for extragalactic sources. 



6. Gamma ray spectrum 

The AW-collisions from the reionization at Zi ~ 5 until today 
give an integrated, somewhat redshifted, gamma spectrum for a 
heavy neutrino with a given mass: 

7 =i C(T) ^r N ^k df dT > (19) 

where C(T) is the clumping factor in Fig. 5 and ^ is the photon 
distribution in Fig. 6. Tq = 2.725 K is the temperature of the 
CMB today and T, is the reionization temperature, which we set 
to Ti = 5 ■ T . 

The resulting E 2 I is presented in Fig. 8. When we com- 
pare the calculated heavy neutrino signal with data from EGRET 
(Sreekumar et al. 1998), we see that only neutrino masses around 
M N ~ 100 or 200 GeV would be detectable, and then only as 
a small bump in the data around E y ~ 1 GeV. For interme- 
diary neutrino masses, the signal would exceed the observed 
gamma ray data. In Fig. 9, the peak intensity for the different 
heavy neutrino masses is plotted, as well as EGRET data for 
the corresponding energy with error bars. The data represent the 
observed diffuse emission at high latitudes (\b\ > 10 degrees), 
where first the known point sources were removed and then the 
diffuse emission in our galaxy was subtracted. 

We have also compared the height of the curves, both with 
and without clumping, and the integrated difference is roughly a 
factor of 30. 

7. Discussion and conclusions 

The numerical calculation of the evolution of the heavy neu- 
trino number density indicates that in the mass region 100 < 
M N < 200, the cosmological neutrinos would give a cosmic ray 
signal that exceeds the measurements by the EGRET telescope 
(Sreekumar et al. 1998). Note that the clumping factor for these 
limits is rather conservative. In Ullio et al. (2002), this factor is 
much larger, which would also produce a stronger limit on the 
heavy neutrino mass. 




M N [GeV] 



Fig. 9. Maximum cosmic gamma radiation from photons pro- 
duced in AW-collisions as a function of neutrino mass (in GeV). 
The marked region is excluded since Q. N > Q.dm within. The 
data are taken at the energy corresponding to the maximum in 
Fig. 8. with error bars. 

We can also compare our neutrino density with the results 
from the Kamiokande collaboration (Mori et al. 1992). We scale 
the neutrino signal in their Fig. 2 to Q/v/Q^m, where we use 
h () = 0.71, n m = 0.2678 and £l b = 0.044. This is shown in 
Fig. 10, where we compare our numerical results for the relic 
neutrino density to the observed muon flux in the Kamiokande 
detector. This gives an exclusion region of 80 < M N < 400 GeV. 
Our analytical results, which are comparable to the traditional 
relic neutrino densities, is about a factor two lower, giving an 
exclusion region of 90 < M N < 300 GeV. The model that gives 
these limits (Gould 1987) is rather complicated and not verified 
experimentally, so these results cannot be taken strictly. Note 
also that in the three-year WMAP analysis (Spergel et al. 2007), 
the value of Q.dm depends on which other data the WMAP data 
are combined with. For WMAP+CFHTLS Q.dm can be as high 
as 0.279 and for WMAP+CBI+VSA it can be as low as 0.155. 
The higher of these possibilities would give an exclusion region 
of 85 < Mn < 350 GeV. The lower boundary value would give 
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Fig. 10. Predicted signal from enhanced AW annihilation in 
the earth and the sun compared to the measured signal in 
the Kamiokande. On the y-axis: the number of muons (per 
100 m 2 year) produced by muon neutrinos resulting from heavy 
neutrino collisions in the sun and the earth, as evaluated by Mori 
et al. (1992), but scaled to our Q,n(Mn). On the x-axis: the heavy 
neutrino mass in GeV. 



an exclusion region of 75 < < 500 GeV. A conservative 
limit based on the Kamiokande data gives the exclusion region 
100 < M N < 200 GeV. 

If a heavy neutrino exists with a mass ~ 100 GeV or 
Mn ~ 200 GeV it would give a small bump in the data at E y ~ 1 
GeV. Currently the data points are too far apart and the error bars 
too large to neither exclude nor confirm the eventual existence 
of such a heavy neutrino. Most of this part of the gamma ray 
spectrum is usually attributed to blazars, which have the right 
spectral index, ~ 2 (Mukherjee et al. 1997). 

We note that there could be an enhancement in the signal due 
to the higher DM densities within galaxies compared to the mean 
density in the halos. On the other hand, from within galaxies 
there will also be an attenuation due to neutral hydrogen, thus 
reducing the enhancement. There will also be a certain degree of 
extinction of the signal due to neutral hydrogen along the line of 
sight, but even if we assume complete extinction above z = 4 the 
resulting spectrum decreases with only about 20%. 

We are also aware of the ongoing debate concerning the an- 
tiprotons - whether or not the DM interpretation of the EGRET 
gamma excess is compatible with antiproton measurements 
(Bergstrom et al. 2006; de Boer et al. 2006). We note the argu- 
ment by de Boer that antiprotons are sensitive to electromagnetic 
fields, and hence their flux need not be directly related to that of 
the photons, even if they too were produced by AW annihilation. 

In the advent of the Large Hadron Collider, we also point out 
that there may be a possibility to detect the existence of a heavy 
neutrino indirectly through the invisible Higgs boson decay into 
heavy neutrinos (Belotsky et al. 2003). 

It will of course be interesting to see the results of the gamma 
ray large area space telescope (GLAST). It has a field of view 
about twice as wide (more than 2.5 steradians), and sensitiv- 
ity about 50 times that of EGRET at 100 MeV and even more 
at higher energies. Its two-year limit for source detection in an 
all-sky survey is 1.6 x 10~ 9 photons cm~ 2 s _1 (at energies > 
100 MeV). It will be able to locate sources to positional accu- 
racies of 30 arc seconds to 5 arc minutes. The precision of this 
instrument could well be enough to detect a heavy neutrino sig- 
nal in the form of a small bump at E ~ 1 GeV in the gamma 



spectrum, if a heavy neutrino with mass ~ 100 or 200 GeV would 
exist. 

There are also some other possible consequences of heavy 
neutrinos that may be worth investigating. The DM simula- 
tions could be used to estimate the spatial correlations that the 
gamma rays would have and to calculate a power spectrum for 
the heavy neutrinos. This could be interesting at least for masses 
M N ~ 100 GeV and M N ~ 200 GeV. The annihilation of the 
heavy neutrinos could also help to explain the reionization of 
the universe. Another possible interesting application of heavy 
neutrinos would be the large look-back time they provide (Silk 
& Stodolsky 2006), with a decoupling temperature of > 10 13 K 
(Enqvist et al. 1989). 

Acknowledgements. E. E. would like to express his gratitude to Konstantine 
Belotsky, Lars Bergstrom, Michael Bradley, Alexander Dolgov, Kari Enqvist, 
Kimmo Kainulainen and Torbjorn Sjostrand for useful discussions and helpful 
comments and explanations. We are both grateful to the GallCS group, who 
has provided the complete dark matter simulations and finally to the Swedish 
National Graduate School in Space Technology for financial contributions. 



References 

Ando, S., Komatsu, E., Narumoto, T., & Totani, T. 2007, Phys. Rev. D, 75, 
063519 

Bardeen, J. M., Bond, J. R., Kaiser, N, & Szalay, A. S. 1986, ApJ, 304, 15 

Becker, R. H., Fan, X., White, R. L., et al. 2001, AJ, 122, 2850 

Belotsky, K., Fargion, D., Khlopov, M., Konoplich, R., & Shibaev, K. 2003, 

Phys. Rev. D, 68, 054027 
Belotsky, K., Fargion, D., Khlopov, M., & Konoplich, R. V. 2004, ArXiv High 

Energy Physics - Phenomenology e-prints, hep-ph/041 1093 
Berezinsky, V., Dokuchaev, V., & Eroshenko, Y. 2006, Phys. Rev. D, 73, 063504 
Bergstrom, L., Edsjo, J., Gustafsson, M., & Salati, P. 2006, J. Cosmology 

Astropart. Phys., 5, 6 
Blaizot, J., Guiderdoni, B., Devriendt, J. E. G., et al. 2004, MNRAS, 352, 571 
Blaizot, J., Szapudi, I., Colombi, S., et al. 2006, MNRAS, 369, 1009 
Cattaneo, A., Blaizot, J., Devriendt, J., & Guiderdoni, B. 2005, MNRAS, 364, 

407 

Coleman, T. S. & Roos, M. 2003, Phys. Rev. D, 68, 027702 

Colless, M., Dalton, G., Maddox, S., et al. 2001, MNRAS, 328, 1039 

de Boer, W., Gebauer, I., Sander, C, Weber, M., & Zhukov, V. 2006, ArXiv 

Astrophysics e-prints, astro-ph/06 12462 
Diemand, J., Moore, B., & Stadel, J. 2005, Nature, 433, 389 
Dolgov, A. D. 2002, Phys. Rep., 370, 333 
Dolgov, A. D. & Zeldovich, Y. B. 1981, Rev. Mod. Phys., 53, 1 
Eke, V. R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263 
Enqvist, K, Kainulainen, K, & Maalampi, J. 1989, Nuclear Physics B, 317, 647 
Fargion, D., Khlopov, M. Y, Konoplich, R. V., & Mignani, R. 1995, Phys. Rev. 

D, 52, 1828 

Fargion, D., Khlopov, M. Y, Konoplich, R. V., & Mignani, R. 1998, Soviet 
Journal of Experimental and Theoretical Physics Letters, 68, 685 

Fargion, D., Konoplich, R., Grossi, M., & Khlopov, M. 1999, ArXiv 
Astrophysics e-prints, astro-ph/9902327 

Gould, A. 1987, ApJ, 321, 571 

Gunn, J. E. & Peterson, B. A. 1965, ApJ, 142, 1633 

Hannestad, S. 2006a, Journal of Cosmology and Astro-Particle Physics, 1, 1 
Hannestad, S. 2006b, Progress in Particle and Nuclear Physics, 57, 309 
Kainulainen, K. & Olive, K. A. 2003, in Neutrino Mass, ed. G. Altarelli & 
K. Winter, 53-+ 

Kolb, E. & Turner, M. 1990, The Early Universe (Oxford: Westview Press), xlii 
+ 547 

Lavalle, J., Pochon, J., Salati, P., & Taillet, R. 2006, ArXiv Astrophysics e-prints, 
astro-ph/0603796 

Maltoni, M., Novikov, V. A., Okun, L. B., Rozanov, A. N., & Vysotsky, M. I. 

2000, Physics Letters B, 476, 107 
Mather, J. C, Fixsen, D. J., Shafer, R. A., Mosier, C, & Wilkinson, D. T. 1999, 

ApJ, 512, 511 

Mori, M., Hikasa, K, Nojiri, M. M., et al. 1992, Physics Letters B, 289, 463 
Mukherjee, R., Bertsch, D. L., Bloom, S. D., et al. 1997, ApJ, 490, 1 16 
Ninin, S. 1999, PhD thesis: Universite Paris 1 1 

Novikov, V. A., Okun, L. B., Rozanov, A. N, & Vysotsky, M. I. 2002, Physics 
Letters B, 529, 111 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, 
Numerical Recipes in C: The Art of Scientific Computing (New York, NY, 
USA: Cambridge University Press) 



8 



Erik Elfgren and Sverker Fredriksson: Mass limits for heavy neutrinos 



Silk, J. & Stodolsky, L. 2006, Physics Letters B, 639, 14 

Sjostrand, T., Mrenna, S., & Skands, P. 2006, Journal of High Energy Physics, 
5, 26 

Spergel, D. N., Bean, R., Dore, O., et al. 2007, ApJS, 170, 377 
Spergel, D. N., Verde, L., Peiris, H. V., et al. 2003, ApJS, 148, 175 
Sreekumar, P., Bertsch, D. L., Dingus, B. L., et al. 1998, ApJ, 494, 523 
Szapudi, I., Bond, J. R., Colombi, S., et al. 2001, in Mining the Sky: Proc. of the 

MPA/ESO/MPE Workshop, ed. A. Banday, S. Zaroubi, & M. Bartelmann, 

249 

The LEP Collaborations. 2001, ArXiv High Energy Physics - Experiment e- 

prints, hep-ex/01 12021 
Ullio, P., Bergstrom, L., Edsjd, J., & Lacey, C. 2002, Phys. Rev. D, 66, 123502 
Yao, W.-M., Amsler, C, Asner, D., et al. 2006, Journal of Physics G Nuclear 

Physics, 33, 1 



